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We show that NS's with large toroidal B- fields tend naturally to evolve into potent gravitational- 
wave (gw) emitters. The toroidal field Bt tends to distort the NS into a prolate shape, and 
this magnetic distortion dominates over the oblateness "frozen into" the NS crust for B t ;> 
3.4 x 10 12 G(i^ s /300 Hz) 2 . An elastic NS with frozen-in B-field of this magnitude is clearly secu- 
larly unstable: the wobble angle between the NS's angular momentum J* and the star's magnetic 
axis re* B grow on a dissipation timescale until J 1 and re* B are orthogonal. This final orientation is 
clearly the optimal one for gw emission. The basic cause of the instability is quite general, so we 
conjecture that the same final state is reached for a realistic NS, with superfluid core. Assuming 
this, we show that for LMXB's with B ; ~2x 10 12 - 2 x 10 14 G, the spindown from gw's is sufficient 
to balance the accretion torque — supporting a suggestion by Bildsten. The spindown rates of most 
millisecond pulsars can also be attributed to gw emission sourced by toroidal B-fields, and both 
these sources could be observed by LIGO II. While the first-year spindown of a newborn NS is most 
likely dominated by electromagnetic processes, reasonable values of Bt and the (external) dipolar 
field Bd can lead to detectable levels of gw emission, for a newborn NS in our own galaxy. 
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I. INTRODUCTION 

Nearly periodic gravitational waves (gw's) from dis- 
torted or wobbling neutron stars (NS's) are among the 
most promising sources for the large-scale gravitational- 
wave detectors now coming online. In this paper we show 
that NS's with large toroidal B-fields tend naturally to 
evolve into potent gw emitters. 

While most known, young pulsars have dipolar mag- 
netic field strengths Bd ~ 10 12 G, studies of soft-gamma 
repeaters and anomalous X-ray pulsars suggest that a 
substantial fraction of NS's are born with larger mag- 
netic fields, B ~ 10 14 - 10 15 G @-§j. In cases where the 
B-field is initially "wound up" by differential rotation (as 
opposed to convective motions), one would expect the 
field's toroidal piece B t to dominate the poloidal part 
Bp. 

The toroidal field B t tends to distort a NS into a pro- 
late shape, and we show that this magnetic distortion 
dominates over the oblateness "frozen into" the NS crust 
for B t £ 3.4 x 10 12 (^/300Hz) 2 G, where v s is the NS 
spin frequency. An elastic NS with frozen-in B-field of 
this magnitude is clearly secularly unstable: the wob- 
ble angle between the NS's angular momentum J 1 and 
the star's magnetic axis n' B grow on a viscous timescale 
until J 1 and n B are orthogonal. This final orientation 
is the optimal one for gw emission. The basic cause of 
the instability seems so general that we conjecture that 
the same final state is reached for a realistic NS (with 
crust and superfluid core). In this paper we consider 
NS's where B t > B p , we assume that most of the results 
for elastic NS's can be carried over to the realistic NS 
case, and examine the consequences. We show that for 
LMXB's with B t ~ 2 X 10 12 - 2 x 10 14 G, the spindown 
from gw's is sufficient to balance the accretion torque- 
supporting a suggestion by Bildsten ||. The spindown 
rates of most millisecond pulsars can also be attributed 
to gw emission sourced by toroidal B-fields of compara- 
ble size, B t ~ 10 12 — 10 13 G. For several known millisec- 
ond pulsars, the gravitational waves would be detectable 
by LIGO II. Finally, while the first-year spindown of a 
newborn NS is most likely dominated by electromagnetic 
(em) processes, reasonable values of B t and the (exter- 
nal) dipolar field Bd can still lead to detectable levels of 
gw emission, for a newborn NS in our own galaxy. 

The idea that a large, toroidal magnetic field should 
drive the magnetic axis orthogonal to J 1 is not new; a 
literature search revealed that P. B. Jones had this idea 
28 years ago M. Jones [fz| — p~0t] went on to develop a sce- 
nario in which the dipole axes of young pulsars (which 
he assumed have B t > B p ) are first driven orthogonal 
to J % by dissipative forces, on a timescale of ~ 10 3 yr, 
but then are slowly driven back parallel to J 1 by em 
radiation reaction, for times J> 10 4 yr. That is, Jones 
estimated that at around age 10 3 — 10 4 yr, dissipative 
forces would diminish to the point that em backreaction 
would drive the dynamics. Clearly, Jones was strongly 



guided by the observation that most young pulsars do 
not appear to have dipole axes orthogonal to J 1 (nor can 
they be perfectly aligned with J 1 , as one would expect for 
Bp > B t , since they pulse). Given his strong argument 
that dissipation should drive the magnetic axis orthogo- 
nal to J 1 , and the apparent contradiction by observation, 
there are three likely ways out: 1) assume dissipative ef- 
fects become negligible on a timescale of order the ages of 
known, young pulsars, 2) assume the pulsar emission re- 
gion is not aligned along a principal axis of the internal 
B-field, or 3) assume that the internal B-field is either 
not born with, or quickly loses, any strong directional- 
ity. Jones developed a model around the first option, 
but found little support for it in the data. 

Here we essentially resurrect the Jones argument, but 
concentrate on a different class of objects possessing 
much smaller external dipole fields Bd'- the LMXB's and 
millisecond pulsars. For these objects, it seems clear that 
dissipation must dominate over the em torque in deter- 
mining the alignment of the NS's axes. We are agnostic 
on the question of why the normal (i.e., young, slow) 
pulsars are generally neither aligned nor orthogonal ro- 
tators. It is irrelevant to our argument whether option 1 
or 2 above is correct; option 3 seems the least likely to 
us, but if option 3 is always correct it would render this 
paper irrelevant. 

Of course, it has been noted before that magnetically 
distorted NS's can generate gw's; see, e.g., Bonazzola 
& Gourgouhlon and D. I. Jones Eg] and references 
therein. Also, Melatos & Phinney |13| recently showed 
how the poloidal B-field on accreting NS's might build up 
"mountains" of accreted matter on the magnetic poles, 
leading to gw emission. What is new here is our observa- 
tion that NS's with large toroidal B-fields should natu- 
rally evolve to configurations with large gw emission, and 
our indentification of LMXB's and millisecond pulsars as 
likely "sites" for this phenomenon. 

In §11 we derive the basic energetics and timescales. 
These are applied to LMXB's and millisecond pulsars 
in §111, and to newborn NS's in §IV. For our estimates, 
we shall always use the following fiducial NS parameters: 
M = 1.4M , R = 10 km, and I = 10 45 g-cm 2 . 



II. ENERGETICS AND TIMESCALES 

Let n B be the NS's magnetic axis (so that if we call 
n l B the z— axis, then the toroidal field B t points in the 
</>— direction). Following Pines & Shaham @, we can 
write the NS's inertia tensor I lJ as the sum of four pieces- 
a spherical piece and three quadrupolar distortions-as 
follows: 
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Here e 13 is the flat, spatial 3-metric, I^e 13 is the spher- 
ically symmmetric part of the inertia tensor, and the 
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terms proportional to en, £d, and eb are the quadrupolar 
distortions due to the star's spin, crustal shear stresses, 
and the magnetic field, respectively. 

The centrifugal piece "follows" the instantaneous 
(unit-)spin vector Uq, but the unit-vectors n d and n l B 
are assumed fixed in the body frame. (However on long 
timescales, n l d can be "re-set" by crustal relaxation.) 

The centrifugal piece en in Eq. is approximately 
given by en w 0.3(^ s /kHz) 2 , where v s is the spin fre- 
quency. The term e^ is due to the elastic crust's "mem- 
ory" of some preferred shape. Absent the magnetic field, 
it is the residual oblateness the NS would have if it were 
spun down to zero frequency, without the crust breaking 
or otherwise relaxing. (To understand this, consider a 
rotating NS with relaxed crust. Removing the centrifu- 
gal force would decrease the NS's oblateness, but then 
shear stresses would also build up that tended to push 
the crust back to its relaxed, oblate shape.) For a NS 
with relaxed crust, e^ is proportional to the centrifugal 
oblateness: td = beQ. The coefficient b was recently cal- 
culated by Cutler et al. fL5| , who solved the coupled hy- 
droelastic equations describing the deformed crust and 
found b w 2 x 10~ 7 . (An older, back-of-the-envelope es- 
timate of b ~ 10~ 5 used by many authors turned out to 
be too high by a factor ~ 40.) Therefore we adopt the 
estimate 
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Again, this estimate assumes the crust is nearly relaxed. 
While this might not be a good assumption for young, 
slow NS's, it seems safe for the fast, old NS's of inter- 
est to us here, for two reasons. First, during the LMXB 
phase, the entire crust is replaced by accreted matter as 
the NS spins up, and it seems unlikely the new crust 
would "remember" the preferred shape from an earlier 
era. Second, the alternative requires the crust to with- 
stand average strains of order 3 x 10 _3 (^ s /300 Hz) 2 for 
~ 10 7 — 10 10 yr without relaxing which also seems 
unlikely. 

Next we need es- Of course, the actual distortion pro- 
duced by some toroidal field B t depends on the precise 
distribution of B-field (and hence currents) within the 
star, which is highly uncertain. So we must be satis- 
fied with a rough estimate of e^, which we derive as 
follows. Let n B be along the z— axis. For an incom- 
pressible, constant-/? star that is not superconducting, 
€b = ijzz — Ixx)/Izz depends only on the total energy in 
the toroidal field 0SM: 



(2.3) 



where Eq = (3/5)GM 2 / R is the star's gravitational 
binding energy. The sign of es is negative because the 
toroidal magnetic field lines are like a rubber belt, pulling 
in the star's waist at the magnetic equator. 



Now, it is generally believed that for interior temper- 
atures T < 10 9 - 10 10 °K, the protons in the NS interior 
are a type II superconductor; hence the magnetic field is 
confined to flux tubes with field strength B c \ w 10 15 G. 
(For B > B c i, the magnetic fluxoids overlap.) Then 
for B < B c \ 1 the anisotropic part of the mean magnetic 
stress tensor is increased over its non-superconducting 
value by a factor B c i/B, and cb increases by the same 
factor (^JTT^l . (The virial-theorem-based derivation of e# 
in Ostrikcr & Gunn |IJ| makes this completely clear.) 
Putting this factor together with Eq. fl2.3|), but using 
E G = (3/4:)GM 2 /R (i.e., the n = 1 polytrope result, 
which quite accurately gives the binding energy of real- 
istic NS's), we adopt the estimate 



-1.6 x 1(T 6 (< B t > /10 15 G) 
-1.6 x 1(T 6 (< B\ > /10 30 G) 



B t < B c i, 
Bt > B cl , 



(2.4) 



where < ... > means "volume-average d ov er the NS in- 
terior." Of course, the first line of Eq. (2.4) only applies 
to NS's older than ~ 0.1 — 1 yr (i.e., old and cold enough 
to be superconducting). 

For simplicity we are assuming, here and below, that 
Bt is sufficiently greater than B p that we can neglect the 
latter. When is this a good approximation? For an in- 
compressible, constant-p, non-superconducting NS, with 
B l of the form Bp^ + Btcj) 1 in the interior and matched to 
a dipolar field in the exterior, the deformation is 1 17 l§fl f\ 
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For a superconducting NS, with B t , B p « B c i 7 we claim 
this becomes 

25i? 4 / \ 
eB = 24GM 2 " [ 2 < BpBcl > ' < BtBcl > ) ■ (2 - 6) 

(The small change (21/10) — > 2 reflects the fact that the 
external poloidal magnetic field energy becomes a negligi- 
ble fraction of the total field energy in the superconduct- 
ing case.) Clearly then, as long as < B t > ^ 4 < B p >, 
our estimates should be fairly reliable. 

We next consider the kinetic energy KE of the spinning 
NS as the orientation of the body changes, for fixed angu- 
lar momentum J 1 . While in the inertial frame the NS's 
angular velocity £l l clearly must remain approximately 
fixed (near the star's angular momentum J*), fi* will mi- 
grate with respect to the body axes as the star precesses 



1 Ostriker & Gunn ac tua lly give 3, not 2.1, as the coef- 
ficient of < Bp > in Eq. (2.6), but we claim they made an 



algebra error; the coefficient on the rhs of their Eq. (A19) 
should be —7/40, not their —1/4. With this correction, their 
coefficent 3 becomes 2.1, in agreement with Wentzel Jl7|. 
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(and as the precession angle secularly changes). Using 
J 1 = P^flj, KE — ^J l fli, and expanding everything to 
lowest order in en, £d and e#, we find: 



KE = 



2I 



1 - en - td{n l d n Ji) 2 ~ ^B(n l B njiY 



(2.7) 



where rij is the unit vector along J\ 

Now consider the case where t B < and \eg\ >> e<2, 
and let 9 be the "wobble angle" between J 1 and n z B . Then 
the ^-dependent piece of KE is (2I Q )^ 1 J 2 (-e B )cos 2 9, 
which is clearly minimized for 9 = tt/2. (More precisely- 
including the effect of non-zero e<j- KE is minimized for 
\n B nji\ = i|(e £ i/es) sm 2/3|, where (3 is the angle between 
n B and n l d .) In addition to the kinetic energy, the wob- 
bling NS stores potential energy (PE) in the stressed 
crust and the distorted fluid. One easily convinces one- 
self that the ^-dependent pieces of PE and KE have 
the same sign and same functional form (oc cos 2 #), and 
that they have comparable magnitudes for NS's with re- 
laxed crusts. Therefore, considering the total energy E 
= KE + PE leads to the same conclusion about energy- 
minimization/stability as considering KE alone. 

It is worth emphasizing that the NS need not actually 
be prolate to be unstable; indeed, for the cases of interest 
here {v s > 100 Hz), the centrifugal oblateness will be or- 
ders of magnitude larger than |e_e | ■ But what matters, for 
both precessional dynamics and stability, is the net pro- 
lateness/oblateness that is "frozen into" the star's body 
frame, not the centrifugal piece that "follows" the star's 
angular velocity. Therefore the condition for instability 
is -e B > e d , or 



B t > 3.4 x 10 12 G 



300 Hz 



(2.8) 



Assuming inequality (2.8) is satisfied, how does the 
instability grow? For an elastic star with frozen-in mag- 
netic field, this is completely clear. In the body frame, 
fl l precesses around n B with precession frequency v s t B 
(neglecting a small correction due to nonzero td)- Due 
to dissipation, the angle between il 4 and n B steadily in- 
creases until they are nearly orthogonal. (For a review 
of NS precessional dynamics, see Cutler & Jones 16 and 
references therein.) 

For a realistic NS, with superfluid core, the preces- 
sional dynamics could be more complicated. But as long 
as those dynamics do not reduce — e B below td (i.e, by 
destroying the large-scale coherence of the internal field 
or building up an equally large internal poloidal piece), 
it seems inevitable, by simple energetics, that dissipation 
will tend to drive n l B orthogonal to J\ In what follows, 
we shall assume this is true. Some support for this as- 
sumption comes from the work of Goldreich & Reiseneg- 
ger |20j| , who show that the composition of a NS interior is 
stably stratified, impeding convective motions that might 
otherwise distort the interior B-ficld. 

We next turn to timescales. There are four relevant 
ones: tqw, tem-, t acc an d 773/5. Once n B has been 



re-aligned perpendicular to J 1 , the NS spins down due to 
gw emission on a timescale tqw = J I {^Jgw), given by 



1/t gw = 5.50 x 10 _1 V 



eg 

10-' 



2 ( v \ 4 
l^kHzJ 



(2.9) 



We must compare this to the timescale tem = 
J I {2 Jem) on which the NS spins down due to electro- 
magnetic processes. Define Bd = 2A4/R 3 , where M. is 
the NS's magnetic dipole moment. (With this definition, 
Bd is the value of the field at the magnetic pole, on the 
NS surface, for an external field that is perfectly dipolar.) 
Then tem is given by 



I /tem 



x 10~ 16 s 



For a ma gnetic dipole rotating in vacuum, the rhs of 
Eq. (2.1C) would be multiplied by a factor sin 2 a, where a 
is the angle between the spin vector and magnetic dipole 
direction. But Goldreich & Julian J2IJ showed that a 
reali stic N S is surrounded by plasma, and as a result 
Eq. ( 2.10| ) is a reasonable estimate of the spindown rate 
even for an ali gne d (i.e., a — 0) rotator. 

From Eqs. fl2.9| )-( 2~T(i| ), spindown due to gw emission 
dominates em spindown for 

B t jB d > 6.2 x 10 3 (^^j , B t <B cl (2.11a) 
./SOOHz^VlO 1 ^ 1/2 



BtjBd > 2.5 x 10 2 



B,, 



B t >B cl (2.11b) 



Such a large ratio B t / Bd would not be expected in a nor- 
mal pulsar, but could occur in a LMXB/recycled pulsar. 
For these sources, accretion apparently "buries" or oth- 
erwise reduces the external dipolar field, but the toroidal 
field in the deep interior could be left largely intact. 

The em torque from the external dipole field causes 
the wobble angle 9 between J 1 and n B to damp (or grow, 
for sin 2 £ > 2/3) on roughly the em spindown timescale 
tem |2 1 : 



1 dsa\9 1 o , 3 9 , , 

-^^r = -2 cos ^-2 sin ^' ^ 

where £ is the angle between the external dipole direction 
and the magnetic axis n l B of the internal B-field. Simi- 
larly, gravitational radiation reaction tends to align the 
distortion axis n' B with the spin direction (independent 
of the sign of e B ) on the gw spindown timescale tqw FH : 
9 = —t^9/t gw for small wobble angle 9, where here t gw 
is defined as the rhs of Eq. ( |2.9| ) . 

Next there is the timescale tacc on which accre- 
tion can significantly change the NS's angular momen- 
tum (either magnitude or direction). We approximate 
tacc = J/2Jacc by H 



1/tacc = 9.3 x lO^V 1 



M 



lO" 9 M /yr 



300 Hz 



(2.13) 
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Lastly, we consider the dissipation timescale tdis on 
which the instability acts. Define n to be the dissipation 
timescale t B is divided by the wobble period, so tjjis = 
nP/e B , where t dis ~ \e B \[J 2 9 2 /(2I)]/E DIS . We can 
rewrite this as 



1/tdis = 3.0 x 10~ 8 s 



10 4 \ / 



V 300 Hz J \10 



(2.14) 



The n factor is hard to estimate, but fortunately most 
of our conclusions are essentially independent of n over a 
huge range. Previous authors who have considered pre- 
cessional damping in NS's have generally concluded that 
precession damps quickly. E.g., Chau & Henriksen |24|] 
considered dissipation in the elastic crust as fi* (and 
hence the centrifugal bulge) precess through the body 
frame. The crust is periodically distorted by the wob- 
bling motion, and the elastic energy stored in the crust 
is some fraction F of the total wobble energy, where we 
estimate 0.05 < F < 0.5. (For the Earth, F is 0.11 gj.) 
If some fraction 1/Q of the elastic energy is dissipated in 
each wobble period, then n ~ Q/F. Typically Q ~ 10 4 
for terrestrial metals, so one might estimate n <; 10 5 . 
Other mechanisms seem likely to lead to faster damping: 
the core and crust will not wobble together rigidly, and 
relative internal motions will lead to frictional damping, 
e.g., at the crust-core interface, or through interactions 
of neutron superfluid vortices with electrons and nuclei 
in the inner crust. (E.g., Alpar & Sauls |26| estimated 
n w 10 2 — 10 4 for slowly rotating NS's, with the dissi- 
pation due to electron scattering off pinned superfluid 
vortices. However superfluid vortices certainly will not 
remain pinned in the wobble regimes of interest to us- 
large wobble angle and v s ;> 300 Hz- so this particular 
estimate does not seem immediately relevant to our case; 
see Link & Cutler 2002 Q.) We therefore leave n as 
unknown parameter to be determined by observation or 
improved theory. The point is that, for LMXB's and 
millisecond pulsars, any n less than ~ 10 9 means the dis- 
sipation timescale is much shorter than all the others, 
and so should successfully drive n B orthogonal to J 1 . 

We have collected all the necessary formulae. We now 
apply them to interesting cases: LMXB's, millisecond 
pulsars, and newborn NS's. 



III. LMXB'S AND MILLISECOND PULSARS 



companion. Accretion reduces the exterior field be- 
low ~ 10 9 G; however in the interior, B t remains in the 
range ~ 10 12 — 10 15 G. As the NS spins back up and its 
external B-field decays, it reaches a state where tjjis is 
shorter than the other timescales (including te m , thanks 
to the decay of B^), while \es\ > (d- Then dissipation 
rapidly "flips" n B perpendicular to J\ Initially n d flips 
along with n B , but we expect the crust of a rapidly spin- 
ning NS will crack/relax such that its new preferred axis 
n d is again aligned with fi*, but now is perpendicular 
to n B . This final configuration is clearly a minimum of 
both KE and PE (for fixed J), independent of the relative 
sizes of ed and \cb\- Thus while the NS continues to spin 
up, dissipation will ensure that n\ remains aligned with 
J 1 , with n B perpendicular to both. This is important 
because td increases quadratically with spin frequency, 
and so (as we shall see) eventually surpasses |es| for the 
fastest LMXB's and millisecond pulsars. Our point is 
that as long as n d has already been "re-set" along J 1 be- 
fore that happens, then n B will remain orthogonal to J\ 
After accretion stops, the NS continues to spin down by 
gw emission. And because n l B is time-varying (rotating 
around J 1 with frequency u s ) while n d is not, it is the 
magnetic distortion that sources the gravitational waves. 

A potential objection to this scenario arises if one as- 
sumes that the pulsar beam is aligned with the axis n l B 
of the B-field in deep interior, since then we would pre- 
dict that (a large fraction of) millisecond pulsars should 
be orthogonal rotators, which is not observed. So for our 
picture to be sensible, we should assume that the B-field 
near the NS surface evolves (e.g., via accretion and crust- 
cracking) in such a way that that the pulsar beam is not 
aligned with n l B . 



B. LMXB's 

The distribution of LMXB spin periods suggest that 
many have "hit a wall" at v s ~ 260 — 590 Hz-far below 
the likely maximum rotation rate of 1 — 1.5 kHz [p|,p8|. 
Bildsten Q proposed that they had reached an equi- 
librium where accretion torque was balanced by spin- 
down from gravitational radiation (resuscitating an idea 
of Wagoner's @). 

Equilibrium between accretion torque and gw spin- 
down implies an ellipticity 



A. Evolution Scenario 



e = 4.5 x 10" 



M 



10- 9 M Q /yr 



1/2 / 300 Hz \ 5/2 



(3.1) 



The arguments and timescales above lead to the fol- 
lowing evolution scenario for (some large fraction of) 
LMXB's and millisecond pulsars. The NS is born with 
B t > 4B p and B d - 10 12 - 10 14 G. At birth, n B is 
likely nearly aligned with the NS's spin (and so is nearly 
aligned with n l d ). The NS spins down electromagneti- 
cally, but much later it is recycled by accretion from a 



For the LMXB's with known spin rates, the range of M's 
is ~ 10~ n — 2 x lO~ 8 M0/yr, with implied e's in the range 
- 3 x 10~ 9 - 3 x 10~ 7 (28). 

On this assumption, gw's from the brightest LMXB's 
(especially Sco X-l) should be detectable by LIGO II. 
Ushomirsky et al. [g3| showed that crustal "mountains" 
can provide the required ellipticity, but only if the crustal 
breaking strain cr max is larger than ~ 10~ 3 — 10~ 2 . 
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Toroidal B-fields place no such requirement on the crust's 
breaking strain: B f ~2x 10 12 — 2 x 10 14 G provides the 
required ellipticity. Our picture requires that gw torque 
domi nates em torque, or Bd < 5 x 10 8 — 5 x 10 10 G, by 
Eq. (2.11a). There are no direct measurements of Bd in 
LMXB's to compare this to, but Bd <; 10 9 G is indicated 
by the lack of X-ray pulses. 

While for simplicity we have assumed B t >> B p , it 
should be clear that the basic picture remains the same 
if they are comparable, so long as the total cb is negative 
and |e B | ~ 10~ 8 - 10~ 7 . 



J1744-1134 (f gw = 491 Hz, r c = 0.72 x 10 10 yr, D = 0.36 



kpc 



yields S/N = 3.6; PSR J1024-0719 (f gw = 387 



Hz, r c = 2.7 x 10 10 yr, D = 0.20 kpc J3j) yields S/N = 
3.3; and PSR J1012+5307 (f gw = 381 Hz, r c = 0.6 x 
10 10 yr, D = 0.52 kpc ||) gives S/N = 2.7. Factor 
— 2.5 improvements in these SNR's could be obtained by 
operating the detector in narrow-band, signal-recycling 
mode. 

Note that since the gw signal is almost completely 
known from radio observations (i.e., up to polarization, 
amplitude and overall phase), S/N J> 2.7 is enough for 
confident detection. 



C. Millisecond Pulsars 

We turn now to millisecond radio pulsars (which we 
define as those with spin periods P < 10 ms), which are 
commonly assumed to be the descendants of LMXB's. 
Our picture suggests than many should now be spinning 
down by gw emission. The NS's characteristic age r e = 
P/(2P) is then just tqw- We can re- write Eq. (|2.9D as 



\e B \ = 6.0 x 10 



-9/ !0 1O yr 



P 



5 msec 



(3.2) 



or B t — 4 x 10 G for our fiducial millisecond pulsar 
values (r c = 10 10 yr and P = 5 msec)-consistent with 
the lower end of the range we inferred for LMXB's. 
(A comparison of derived e^'s in the LMXB and mil- 
lisecond pulsar populations would be interesting, but is 
clearly complicated by selection effects; e.g., low-|e_B| 
millisecond pulsars live longer, and so will dominate 
the observed sample.) Note that \cb\ > £d only for 
P £ 4.0msec(T c /10 10 yr) 1 / 2 (from Eqs. fO and [O), so 



the oldest/fastest millisecond pulsars would not be unsta- 
ble in the way we described in §11. However, as explained 
in §111. A, that probably does not matter: their magnetic 
axes were flipped orthogonal to J 1 while they were being 
spun up, and will not flip back if the crust has relaxed. 

Given a millisecond pulsar at distance D, spinning 
down by gw emission, it is straightforward to compute 
the dctcctability of the gw signal h(t). The rms S/N 
(averaged over source direction and polarization) is 



S/N 



i.i(%£ 



1/2 /lfjiOyr V /2 (2 x 10~ 24 



(3.3) 



where / = 2v s is the gravity- wave frequency, Sh(f) is 
the (single-sided) noise spectral density, and To is the 
observation time. The broad-band LIGO-II noise curve 

1/2 

has a broad minimum at f gw — 400 Hz, with S h ' (/) f» 

2 x 10~ 24 there [0. Using this fiducial value for S^ /2 (/), 
ones finds there are at least 4 known millisecond pulsars 
with S/N > 2.7 in a 1-yr observation. PSR J0437-4715 
(f gw = 347 Hz, t c = 0.6 x 10 10 yr, and D = 0.18 kpc §§) 
gives the highest S/N ratio: S/N = 7.9. Also, PSR 



IV. NEWBORN PULSARS 

As another application of our basic idea, consider the 
spindown gw's from a supernova in our own galaxy (D — 
10 kpc). (The supernova rate in our galaxy is — l/40yr, 
so this author expects to be around for the next one.) It 
does not seem outlandish to posit field strengths Bd — 
10 14 G and B t — 10 15 G. Now em radiation dominates the 
spindown torque, but there is significant gw emission as 
well. 

The matched- filtering (S/N) 2 for a NS spinning down 
solely due I — m — 2 gw's is [B0J 



S/N =- 



2G 



rD 2 c 3 



fma * dJ/df 



fS h (f) 



df • 



(4.1) 



For our case, we just need to multiply the numerator in 
Eq. (LI) by Jgw/Jem = t E m/t gw oc f 2 . Approximat- 
ing the high-frequency part of the broad-band, LIGO-II 
noise spectrum by S h (f) = 2.1 x 10- 47 s(//10 3 Hz) 2 |l), 
we find 



D 



10" 



V B d 



fn 



| 1/2 



(4.2) 



Here we can take f m in to be the frequency (~ 500 Hz) 
where the actual LIGO-II broad-band noise curve "flat- 
tens out." 

E.g., consider a NS with B t = 2 x 10 15 G and B d = 
10 14 G, which spins down from v s — 1 kHz to v s — 250 
Hz in — 10 days (and to v s w 80 Hz after one year). 
In order for n B to "flip" before the NS spins down sub- 
stantially, we needrp/s < tem, or n <J 1.3 x 10 3 , by 
Eqs. (|2. lOf ) and (|2.14|) . Assuming n is this small, and 
taking [Hf max /f min )}^ 2 « 1, we find S/N = 75 for 
D = 10 kpc, which would likely be detectable even in a 
semi-blind search-i.e., one where neutrino detectors and 
other channels give the collapse time and a rough posi- 
tion, but where there is no radio signal to give us the 
rotational phase as function of time. 
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